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ABSTRACT 

The search for planets around White Dwarf stars, and evidence for dynamical instability around them in the 
form of atmospheric pollution and circumstellar discs, raises questions about the nature of planetary systems 
that can survive the vicissitudes of the Asymptotic Giant Branch (AGB). We study the competing effects, on 
planets at several AU from the star, of strong tidal forces arising from the star's large convective envelope, and 
of the planets' orbital expansion due to stellar mass loss. We, for the first time, study the evolution of planets 
while following each thermal pulse on the AGB. For Jovian planets, tidal forces are strong, and can pull into 
the envelope planets initially at ^ 3 AU for a 1 Mq star and ^ 5 AU for a 5 Mq star Lower-mass planets 
feel weaker tidal forces, and Terrestrial planets initially within 1.5 — 3 AU enter the stellar envelope. Thus, 
low-mass planets that begin inside the maximum stellar radius can survive, as their orbits expand due to mass 
loss. The inclusion of a moderate planetary eccentricity slightly strengthens the tidal forces experienced by 
Jovian planets. Eccentric Terrestrial planets are more at risk, since their eccentricity does not decay and their 
small pericentre takes them inside the stellar envelope. We also find the closest radii at which planets will be 
found around White Dwarfs, assuming that any planet entering the stellar envelope is destroyed. Planets are in 
that case unlikely to be found inside ^ 1.5 AU of a White Dwarf with a 1 M0 progenitor and ^ 10 AU of a 
White Dwarf with a 5 Mq progenitor 

Keywords: planets and satellites: dynamical evolution and stability — stars: AGB and post-AGB — stars: 
evolution 



1. INTRODUCTION 

The discovery of White Dwarf stars with metal-polluted 
atmospheres or with circumstellar discs strongly suggests 
that planets, feeding asteroidal or cometary material to lo- 
cations close to the star, exist around these de generate stel- 
lar remnan ts (Zuckerman et al. 2003; Gansicke et al.l 120061 : 
iFarihi et arii2009.: Gansicke et al..,2012) . Furthermore, there 
are claims of planetary detection, by means of variations in the 
timing of stellar pulsations or binary eclipses, around various 
kinds of evolved star (e.g., MuUallv et al. 2008; Oian et al. 
|2009). For such planets to remain to late times they must, 
however, first survive the late stages of stellar evolution. Dur- 
ing the giant phases of the host star, the star loses a large frac- 
tion of its mass, causing the planet's orbit to expand in order 
that angular momentu m be conserved ( Hadiidemetrio ull 1 9631: 
Alexander et al."1976; Livio & Sokenil984t 1 Villaver & Livid 
2007, 2009; Veras et al. 2011). Countering this, however, the 
stellar envelope expands to a radius of several AU, and tidal 
forces wax stronger, potentially drawing plan ets in to their de- 
struction in the envelope (Rasio et al. 1996; Villaver & Livid 
120071 l2009i) . Dynamically, the survival of planets with semi- 
major axes of a few AU is governed by the balance between 
the tidal force which pulls the planet towards the expanded 
envelope, and the effects of stellar mass loss which push the 
planet away. 

The fate of planets around Red Giant Branch (RGB) 
stars has received much attention, since many planets have 
been discovered around sub-giant, RGB and horizontal 
branch stars (e.g.,^ iedzielski et al. 2009; Bowler et al. 2010; 
iJohnson et al.l 1201 It Wittenmver et al. 2011; a recent com- 
pendium is given in fCettel et aL,2012) . Such stars afford the 
only way to probe, with radial velocity surveys, the popula- 
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tion of planets at radii of up to a few AU around stars more 
massive than the Sun, and a lack of close-in planets around 
these stars has been found. The cause of this deficiency has 
been attributed bot h to differen t planet formation processes 
around these stars (ICurriell2009l) . and to tidal engulfment of 
close-in planets as the st ar expands and ascends the RGB 
jSato e t al. 2008?, iVillaver & Livio 2009). Regardless of the 
planetary population on entering the RGB, studies find that 
after the RGB phase planets around 1 Mq stars should have 
survived if they are located beyond a few AU, while around 
more massive stars the minimum radii at which planets should 
have survived are smaller, since the more massive stars have 
much smaller maximum radii on the RGB than Solar mass 
stars (Sato et al. 2008; Villaver & Livio 2009; Kunito mo et alj 
1201 lb . The survival radius increases with pl anet mass as tidal 



force s are stronger for more massive planets (I Villaver & Livid 
120091) . There should, therefore, be a large population of plan- 
ets at radii of several AU which will survive until the asymp- 
totic giant branch (AGB) phase of evolution. A few plan- 
ets, moreover, have bee n found in tight orbits around stars 
on th e horizontal branch ([Silvotti et al]|2007t ICharpinet et alj 
120111: ISetiawan et al.l 120 lOl) . These planets could have sur- 
vived the evolution of the star along the RGB, in the process 
possibly trigg ering mass lo ss to f orm an sdB star dGeier et alj 
l2009l lHebdl2009t [Soke] [HH), or could be the remnants 
of the tidally destro yed metallic cores of massive planets 
(iBear & SokerilMl) . 

The population of planets around White Dwarf stars is, 
however, completely unknown. Despite the advantages that 
White Dwarfs offer to direct detection of massive planets in 
the infrared (Burleigh et al. 2002), and the heroic efforts in 
the search for these plane ts using this technique (Hogan et al] 
^009t ISteele et alj l20Tlh and others CMuUallv et al. 2008 



Debes et alj201 ItlFaedi et al.ll201 Ih . no planet has so far been 
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confidently detected orbiting a White Dwarf. Studies of plan- 
etary dynamics during the preceding AGB phase can therefore 
play a useful role by helping to inform such searches by pre- 
dicting the properties of planets that survive the AGB. 

On the AGB is attained, for many stars, the greatest mass 
loss; this, coupled with the large stellar radius and therefore 
strong tidal forces, means that the AGB phase is likely to be 
the key determinant in the survival of planets to the WD stage. 
The survival of planets during this phase has, however, re- 
ceived less attention. The survival of giant planets during and 
after the AGB was investigated by Villaver & Livio (2007), 
in which study were included the effects of stellar mass-loss, 
the expanding stellar envelope, and evaporation of the planet 
during the Planetary Nebula phase. Tidal effects were, how- 
ever, ne glected. These were included in a more recent study 
(iNordh aus et al. 2010), which followed the evolution of giant 
planets and brown dwarfs from the ZAMS to the WD phase 
under tidal evolution and stellar mass loss using the Reimersi 
(fT975h prescription. 

In the present paper we extend previous work on the fate of 
planets around AGB stars, following the planets' orbital evo- 
lution under the opposing effects of stellar mass loss and tidal 
forces. We make use of realistic stellar evolutionary models 
that follow the star's evolution through each thermal pulse, 
during which the stellar radius experiences large variations, 
and that have a mass loss prescription tested to be valid during 
this phase. This enables us to calculate the changing strength 
of the tides as the star evolves. We extend the previous studies 
to planets of lower mass, which can be dynamically signifi- 
cant on the WD stage in their interactions w ith remnant plan- 
etesimals ( Bonsor. Mustill. & Wvatlil201 Ih . Also considered 
are the effects of including a moderate planetary eccentricity, 
and of changing the several parameters in the tidal formalism. 

This paper is organized as follows. In Section|2]we set out 
the equations governing the evolution of the planet's orbit. In 
Section [3] we describe the stellar models used. The results 
of our numerical integrations are presented in Section |4] and 
discussed in Section|5] We conclude in Section|6] 

2. GOVERNING EQUATIONS 

We adopt the formalism of lZahnI (119771) for the tidal forces 
raised on the star by the planet. The dominant tidal effect 
is the damping of the equilibrium tide by viscous dissipation 
in the star's large convective envelope, which spans almost 
the entire stellar radius ( see Section O. It has been shown 
dVerbunt & Phinne"vll 19951) that this is a good model of tidal 
forces in stellar binaries with a giant component, but it is un- 
clear whether it will be accurate for the weaker tides raised by 
planets. We also consider tidal dissipation in the planet, but 
shall find that it is unimportant. We include the leading-order 
eccentricity contribution to the tidal equations. The contribu- 
tion of the stellar tide to the rates of change of orbital elements 
is: 



(3) 

Here, a and e are the semi-major axis and eccentricity of the 
planet; Af^, Mem and il/pi are the masses of the star, stel- 
lar envelope, and planet; Ri, is the stellar radius; iconv is the 
convective time-scale in the stellar envelope, given by 



tc. 



1/3 
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being the stellar luminosity, i?onv the radius of the stellar 
envelope, and rjp a parameter of order unity. The coefficients 
f i in the tidal equations contain the frequency dependence of 
the components of the tidal forces, and are given by 



fi = /'min 



1, 



27r 



7f 



(5) 



(7i — in, n being the mean motion, are the individual 
frequency components. /' and cp are order-unity parame- 
ters, and 7f an exponent determining the nature of the fre- 
quency dependence. Anal ytical and numerical work suggests 
that -7F sa 1 - 2 CZahn. J977t iGolfreich & Nicholsoni ri977^ 



iGoodman & Ohlll997l: iPenev et al.ll2007l) For compatibihtv 
with the results of IVillaver & Livic (120091) . we adopt rjp = 3, 
cp = 1, /' = 9/2 and 7f = 2 as our fiducial parameters. 

While we have chosen to study the effects of a particular 
tidal mechanism, we recognize that others may be at work. 
These could be dealt with by using the well-known Q formal- 
ism, which has the advantage of putting the unknown strength 
of the tidal forces into the single parameter Q^. However, the 
value that this should take is not certain for giant stars. Stud- 
ies of Main-Se quence pl anet hosts suggest in the range 
lOS-iQio (-e.g.. lJackson et al. 2008; Penev & Sasselov 201 1). 
However, even this poor calibration may not be applicable to 
giant s tars with their very different structure: INordhaus et alj 
(2010) find that their use of the Zahn formalism gives equiva- 
lent Q'^ values of order lO'^ — IQ-^. Use of a Q formalism with 
a Main-Sequence calibration will result in drastically weaker 
tides, while a calibration for giants must come from estimat- 
ing the tidal strength per the assumed dissipation model, with 
little change if the turbulent dissipation model be adopted. 
Other models proposed for objects with c onvective envelopes, 
e.g., d issipation of inertial waves (e.g., Goo dman & Lackne^ 
|2009|), will likely give much higher Q values than the turbu- 
lent dissipation mechanism. 

We also include, for the eccentric planets, a planetary 
tide. The dominant tidal mechanism in giant planet atmo- 
spheres is not known, with dynamical tide mechanisms such 
as gravity wave damping (Lubow, Tout, & Livio 1997) and 
inertial wave damping (Ogilv ie & LinI 120041: IWul l2005aibl : 
IGoodman & Lackner 2009) having been proposed. In light 
of this uncertainty, we adopt the sta ndard Q formalism and a 
constant time lag of the tidal bulge dMatsumura et alj|2010l) . 
This affects the evolution of the orbital elements as 
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with Api = in 1 2(5p[ 1 2ripi — n | being the product of the mean 
motion, the lag time, and the planetary Love number. A con- 
stant time lag entails varying Qpj; our Qpj values quoted be- 
low are normalized to 0.05 AU for comparison with planets 
around more commonly studied Main Sequence stars. We test 
a range of Qpj values, from lO'^ to 10^° for giant planets and 

10" to 10'' for terrestrials; however, we find that in all cases 
the planetary tide is so weak as to have no discernible effect 
on the orbital evolution. The planetary angular momentum 
ilpi also evolves: 



Mi 



pi 



Mpi (M, + Mpi) 



2 ^pl" 



(8) 
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The star is assumed to be non-rotating, which is likely a rea- 
sonable approximation as the spin rates of giants are very low: 
rotational velocities of RGB stars are on average v sin i < 
2km (see e.g. lde Medeiros & Mayol^l999^ . with rapid ro- 
tators representing less than 2% of a sample of 1300 K giants 
(^arlberg et al. 2011). Conservation of angular momentum 
should ensure that stars during the AGB phase are slow ro- 
tators as well; for they will spin down after the horizontal 
branch as their radii expand. 

Finally, the planet's orbit can expand as the star loses mass. 
All expansion is assumed to be adiabatic, with the eccentricity 
remaining constant and the semi-major axis varying by 



'Af, + Mpi' 



(9) 



although we do not include the rate of change of the 
planet's mass due to ei ther mass accretion or evaporation (see 
IVillaver & Livid 1200 9') because it is negligible compared to 
the stellar mass loss term during the AGB phase. The as- 
sumption of adiabaticity is justified so long as the time-scale 
for mass loss is much longer than the Keplerian time -scale. 
This criterion was parameterized bv lVeras et al.l (1201 Ih as 
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27rMoyr-i VAU/ \Mq 



-3/2 



(10) 



finding that substantial deviations from adiabaticity could 
arise for "if as small as 0.02, this causing eccentricity varia- 
tions of order 0. 1 if maintained for several thousands of years. 
For our stellar models, taking a planet at 5 AU, the maximum 

is 3 X 10^'^, attained only briefly, so we expect that the 
adiabatic approximation will hold. 

We negl ect any drag forces f rom th e strong stellar wind, 
shown by Duncan & Lissauerl (119981) to reduce a planet's 
semi-major axis by a mere 1 per cent over the AGB Hfetime. 

3. STELLAR MODELS 

As described previously, the chosen parameterization to 
calculate the tidal forces involves the dissipation of energy 
in the convective envelope. In order to calculate its effect on 
the orbital evolution of the planet a number of parameters re- 
lated to the stellar interior such as the mass and radius of the 
convective envelope are needed. A good parameterization of 
the mass loss during this phase is especially critical as well, 
because its time-scale is linked to the stellar structure and be- 
cause it is a determining factor in the planet's orbital evolu- 
tion. 



The internal structure of an AGB star is complex, espe- 
cially during the thermally-pulsing stage when strong time- 
dependency occurs and when most of the mass loss takes 
place. During this phase two nuclear burning shells, the in- 
ner thermally unstable and burning helium, the outer com- 
posed o£JiYdrogenj_suiTouii^ carbon-oxygen 
core (ISchwarzschild & HarmI [T965I) . During each thermal 
pulse the inner edge of the convective envelope moves in- 
wards, close to the H-burning shell, dredging up fresh ma- 
terial to the surface. Variations in the surface luminosity and 
radius are produced as a consequence of the thermal read- 
Justments taking place dur ing the thermal pulses (see, e.g., 
iKarakas & Lattanzioll2007l) . 

These thermal pulses responsible for the modulations in the 
stellar radius are also the drivers of the mass loss through 
the generation of shocks in the stell ar envelope tha t allow the 
nucleation of dust grains (see, e.g.. lBowenl ll988'. and refer- 
ences therein). As grains nucleate and grow they experience 
the force exerted by the stellar radiation pressure and thus are 
accelerated away from the star. The momentum coupling be- 
tween gas and dust then drives the mass outflow. 

We have used the temporal evolution of mass loss provided 
by the stellar evoluti onary rn odels of Vassiliadis & Woo^ 
d 19931) : see also Vill aver et aH (l2002l 1201 2l) for a discussion 
on the mass-loss in these models. For consistency, the stellar 
structure is derived from the same models using the following 
approximations. First, the envelope mass Mem is approxi- 
mated by Ali, — Mc, where Af* is the star's total mass, and Afc, 
the core mass, has been taken as its value at the first thermal 
pulse from the parameterization, a s a funct i on of the initial 
mass and metallicity, given by Kara kas et al.l(l200l . The core 
mass is expected to increase slightly due to H-burning during 
the inter-pulse periods. This increase in mass, although funda- 
mental in problems devoted to stellar evolution, is negligible 
for our purposes, being a few hundredths of a Solar mass. The 
use of the initial core mass, so that the envelope mass never 
goes to zero, also avoids numerical problems when integrating 
the equations of motion (Equations[T]|2]and|5]l. The envelope 
mass is then made to decrease with the stellar mass. Second, 
we have taken the radius of the convective envelope to be the 
radius of the star. The region interior to the H-burning shell is 
similar in size to a white dwarf, with almost the entire stellar 
volume taken up by the convective envelope, and the radius at 
the base of the convectiv e envelope being just « 10^"* i?^ (see, 
e.g.. lMarigo et alJll998l) . This approximation thus has a neg- 
ligible effect when it enters in the estimate of the convective 
time-scale in Equation (HJi. 

4. NUMERICAL STUDIES 

4.1. Planets on circular orbits 

We first considered the evolution of planets initially on cir- 
cular orbits about the stars of initial mass 1, 1.5, 2, 2.5, 3.5 and 
5M0, for the duration of the stars' thermally pulsing AGB 
lifetime. Note that the "initial" mass is that at the beginning 
of the AGB; this may have been reduced below the MS mass 
if mass loss has already occurred on the RGB. We considered 
planets of Terrestrial (IM©), Neptunian (17.1 M®) and Jo- 
vian (318 M®) mass. The tidal equations in Section|2l includ- 
ing orbital expansion fro m mass loss, were integrated with a 
Runge-Kutta integrator ( Press et alj|2007h . The integration 
was stopped if the planet at any time entered the stellar enve- 
lope. The lifetime of planets having entered the envelope is 
very short, and despite the brevity of some of the radial pulses 
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Figure 1. Orbital evolution of planets initially at different semi-major axes 
around a 1 Mq AGB star, under the effects of tidal forces and stellar mass- 
loss. Planets in the top panel are of Jupiter's mass; those in the bottom are of 
Earth's mass. Black lines show planets that survive; grey lines (red on-line) 
show planets that strike the stellar envelope. 

we expect that the orbits of such planets will decay rapidly. 
Loss of orbital ener gy is very fast and the time- scale can be 
as short as weeks (Nordhaus & Blackman'2006^. The ques- 
tion of the potential final fate of planets engulfed in the stellar 
envelope is not straightforward. IVillaver & Livid (2007) esti- 
mated that a planet up to 15 Mj can be evaporated inside the 
envelope of an AGB star with a main-sequence mass of 1 M0, 
with this mass limit increasing and reaching well into the 
stellar- regime 120 Mj, 0.11 M©) if the planet (or brown 
dwarf) is engulfed in the envelope of a 5 M© AGB star How- 
ever, as has been argued sever al times (Villaver 2011, 2012!) 
the value of this maximum mass is very uncertain because it 
depends on several factors, such as the efficiency of stellar en- 
velope ejection, which are largely unknown. Furthermore, to 
survive the Common Envelope phase, a planet must be able 
to supply enough of its orbital energy to the stellar envelope 



Figure 2. As Figure[T]but for a 1.5 Mq star. Note here, and also in subse- 
quent plots, the different axis scales. 

to unbind the envelope before it spirals down t o the Roche ra- 
dius, at which the planet is tidally disrupted. Nordhau s et alj 
(i2010.) estimated the minimum mass for this to happen at just 
under 10 Jovian masses for a 1 star, much larger than the 
planets we are considering. Note that the Roche radius is very 
small compared to the radius of the star, so planetary disrup- 
tion following passage of the Roche radius will only occur if 
the planet has already entered the envelope. 

In Figures [T] to |6] we show the semi-major axis evolution 
of planets of Jovian and Terrestrial mass around the differ- 
ent stars; the stellar radius is also shown. The evolution of 
the planets orbiting each star displays the expected conflict 
between orbital expansion due to stellar mass loss and or- 
bital contraction due to tidal effects. Towards the end of each 
thermal pulse, tidal decay accelerates as the stellar radius ex- 
pands, but may then be reversed as the star loses mass at the 
end of the pulse. Some planets undergo several cycles of mo- 
tion inwards and outwards before finally colliding with the 
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Figure 3. As Figure[T]but for a 2 Mq star. 

Star or escaping to safe radii. The evolution of the planet's 
orbit depends only on its mass and initial semi-major axis, 
with all planets interior to a critical radius being engulfed and 
all exterior avoiding engulfment. For brevity, we refer to the 
engulfed planet with the largest initial semi-major axis as the 
moriturus ultimus. The orbit of the moriturus ultimus is at 
several AU and is larger for the more massive stars and more 
massive planets. These radii are plotted in Figure [T] as black 
lines. 

For the Jovian planets, tidal decay can be significant, since 
tidal forces are strong for massive planets. Around lower mass 
stars, tidal forces can pull in planets that have initial orbits far 
larger than the maximum stellar radius. This is not the case, 
however, for highest mass star (5 M0), where the different 
pattern of stellar radius evolution and mass loss, and smaller 
planet: star mass ratio, means that some planets which begin 
interior to the maximum stellar radius will be moved beyond 
it by mass loss before this radius is attained, and be thus saved 
from entering the stellar envelope. 



2.5 Mq 




1 .... I .... I .... I .... I .. I . I 
0.0 0.5 1.0 1.5 2.0 2.5 
time [Myr] 



2.5 Mq 




time [Myr] 
Figure 4. As Figure[T]but for a 2.5 M© star. 

For the Terrestrial planets, tidal forces are very weak, and 
the evolution is dominated by the expansion due to mass-loss. 
This is in all cases sufficient to push some planets that begin 
interior to the maximum stellar radius out to safety. This is 
particularly noticeable in for the 5 progenitor, which ex- 
pands to over 5 AU, planets beyond 2.8 AU however being 
safe. The entry of these planets into the stellar envelope typi- 
cally occurs when the stellar radius suddenly swells to engulf 
them. 

Neptune-mass planets, as shown in Figure [T] display be- 
haviour intermediate between Jovian and Terrestrial planets, 
as is to be expected. The semi-major axis of the moriturus 
ultimus for each star is slightly larger than for the Terrestrial 
planet, and in all cases save the 1 star it lies within the 
maximum stellar radius. 

4.2. Inclusion of orbital eccentricity 

We now turn to consider the effects of moderate planetary 
eccentricities on the evolution of the planets' orbits around 
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Figure 5. As Figure[T]but for a 3.5 Mq star. 

AGB stars. Because the ratio of the planetary radius to the or- 
bital radius is so small, the eccentricity evolution is governed 
entirely by the stellar tide, in contrast to most cases of planets 
orbiting Main Sequence stars where both planetary a nd stel lar 
tides are effective (e.g., Jackson et al. 2008; Hansenl l2010l) . a 
distinction pointed out bv jNordhaus et al.i (|201Q|) . Here, we 
include the lowest-order terms in eccentricity in the equations 
for the stellar tides as described in Section|2l and also include 
for completeness the lowest-order contribution to the plane- 
tary tide as described in Section|2] Because of the depen- 
dence of the tidal response on forcing frequency, an expansion 
to arbitrarily high eccentric ity using the Eggleton formalism 
(lEggleton. Kiseleva. & Hu t 1998) is inconsistent with our as- 
sumed tidal model. We therefore restrict attention to planets 
with an initial eccentricity of 0.2, for which the lowest or- 
der equation s were found to be acc urate when compared with 
the complete lEggleton et al.l(ll998h equations for the constant 
time-lag model. 
The tidal evolution of example planets orbiting a star of 
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Figure 6. As Figure [T]but for a 5 Mq star. 

1 M0 is shown in Figure[8] That of Jovian planets is shown in 
the top panel. The eccentricity and semi-major axis damp on 
similar time-scales, and the inclusion of eccentricity shortens 
the lifetime of the Jovian planet by about 10 per cent com- 
pared to the lifetime of a non-eccentric planet. For other 
initial conditions, reductions of over 20 per cent were seen. 
The eccentricity has decayed by a factor of 10 by the time the 
planet hits the stellar envelope. The evolution of planets with 
quality factors Qpj — 10^ and 10^° was considered, com- 
fortably bracketing estimated Q values for gas giant planets, 
but planetary tides are so weak that they are indistinguish- 
able. The bottom panel of Figure [8] shows the evolution of 
Terrestrial planets, with the two Q'^^ values considered, 10'' 

and 10'^, again being indistinguishable. Here, stellar tides too 
are weak: the eccentricity barely decays, the inclusion of the 
eccentricity only slightly enhances the semi-major axis decay, 
and yet the lifetime of the planet is reduced by a third. This 
is because the planet's pericentre takes it inside the envelope 
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Figure 7. Maximum radius of the stai' on tlie AGB (solid line) and the orbital 
radii of the initially most distant planet engulfed by the star, as a function of 
stellar mass. Black lines show planets on circular orbits; grey lines (red in the 
on-line version), planets initially on eccentric orbits. 

during an earlier radial pulse than the one which swallows 
the planet on a circular orbit. For other initial conditions, the 
planet's lifetime could be halved by giving it an eccentricity 
of 0.2. 

We now determine the radii of the morituri ultimi for ec- 
centric planets. In Figure |9] we show the fates of suites of 
eccentric planets around a 5 Mq star. The top panel shows 
the evolution of Jovian planets. These experience significant 
eccentricity decay prior to engulfment. The bottom shows the 
evolution of Terrestrial planets. These experience very little 
eccentricity or semi-major axis decay, and their fate is de- 
cided by their small pericentres bringing them within the en- 
velope. The radius of the moriturus ultimus is larger when the 
effects of eccentricity are included, by up to 0.25 AU for the 
Jovian planets, up to 0.27 AU for the Neptunians, and up to 
0.6 AU for the Terrestrials. In the first case this is attributable 
to the enhanced semi-major axis decay when eccentricity is 
included; in the last, to the pericentres' being significantly 
smaller than the semi-major axis; both factors contribute for 
the Neptunians. All the critical radii, of both circular and ec- 
centric planets, are shown in Figure [T] 

4.3. Planets that avoid engulfment 

It is of great interest for predicting the nature of planetary 
systems around WD stars to determine the final radius of the 
planet closest to the star that evades engulfment, a planet we 
call the elapsus citimus. These radii are shown in Figure [TOl 
They are larger for larger progenitor mass, both because larger 
stars lose a greater fraction of their mass on the AGB, caus- 
ing greater adiabatic orbital expansion, and because the larger 
maximum radii of such stars cause the engulfment of more 
distant planets. For the Jovian planets, tidal effects are strong 
enough to strand some planets so that they finish up inside the 
maximum stellar radius. This ef fect is larger the s maller the 
mass of the star, as seen also in IVillaver & Liviol (12009') for 
the RGB. Such stranding requires a good deal of fine-tuning 
however, and is infrequent. For this reason, it is also possible 
that we missed some closer-in planets with our sampHng of 
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Figure 8. Top: Tidal evolution of Jovian planets ai'ound a 1 Mq AGB star. 
The solid black line shows the stellar radius; the dashed black line the semi- 
major axis of a planet on a circular orbit. The solid grey line shows the semi- 
major axis of planets with initial eccentricity of 0.2. There are two lines here, 
with Qj^j = 10^ and 10^°, but they are indistinguishable. The dotted grey 
hne (red on-line) and right-hand axis show the evolution of the eccentricity. 
Bottom: the same, for Terrestrial planets. Here the two eccentric planets have 
quality factors 1 and 1000, but again their evolution is indistinguishable. 

parameter space (an initial separation of 8 x 10^"* AU). How- 
ever, due to the non-smooth nature of the stellar radius evolu- 
tion, it is not guaranteed that the elapsus citimus can end up 
arbitrarily close to the star. A planetary eccentricity decreases 
the radius of the Jovian elapsus citimus: it is pulled in more 
because of the stronger tidal forces. 

We show the effects of tidal forces on shrinking the orbit of 
the elapsus citimus in Figure [TT] which shows the ratio of the 
final orbit of each surviving simulated planet, of Jovian mass 
and on a circular orbit, to that expected purely from adiabatic 
mass loss. Any planet lying within the dotted line of Figure|2] 
is engulfed by the envelope and not shown. Just outside this 
radius, tidal orbital decay competes with the orbital expansion 
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Figure 9. Semi-major axis evolution of eccentric planets around a 5 Mq star. 
Initial planetary eccentricity is 0.2. Top: Jovian planets. Here the eccentricity 
decays to a small value prior to engulfment. Bottom: Terrestrial planets. Here 
there is little eccentricity decay. 



due to mass loss, and the orbit of the planet is reduced much 
beyond that expected from pure adiabatic expansion. For this 
reason, the radii of the elapsi citimi plotted in Figure [TO] are 
not simply the radii of the morituri ultimi scaled up by a factor 
of the fraction of the stellar mass lost. However, for planets 
with initial orbits beyond around half an AU greater than those 
of the morituri ultimi, tidal forces are weak and adiabatic or- 
bital expansion is a good approximation. 

For the Earth-mass and Neptune-mass planets, tidal effects 
are weak enough that no stranding inside the maximum stellar 
radius occurred. A lower planet mass means that the elapsus 
citimus ends up more distant from the star, although it begins 
closer than for a more massive planet. This is because the 
larger planets experience stronger tidal orbital decay. The ra- 
dius of the Terrestrial elapsus citimus is only just short of the 
initial radius of the moriturus ultimus scaled up by adiabatic 



Figure 10. Final radii of the innermost surviving planet, along with the max- 
imum stellar radius on the AGB. Lines are as in Figure [T] 
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Figure 11. Ratio of planets' final semi-major axis to that expected from adi- 
abatic orbital expansion, as a function of initial semi-major axis. Tidal forces 
cause the final axes of many planets that are not engulfed to be significantly 
below those expected from pure adiabatic orbital expansion. Sets of points 
are in order of increasing stellar mass, from left to right. 

mass-loss. For the low-mass planets, a planetary eccentric- 
ity means that the elapsus citimus is at a greater radius than 
for the circular planets, in contrast to the case for the Jovian 
planets, because the planetary eccentricity does not add sig- 
nificantly to the weak tidal forces. While post- AGB stars with 
1 M0 progenitors will be cleared of planets out to two or three 
AU, those with 5 Mq progenitors will be cleared of Jovian 
planets out to 10 AU and of terrestrials out to 15 or more AU. 

4.4. Sensitivity to parameters 

There are four parameters in our tidal model which are 
poorly known. These are: cp, governing the onset of the 
frequency dependence of the tidal response, and 7f its expo- 
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nent (Equation |5]); ryp, used to estimate the convective time- 
scale (Equation lU; and /' which governs the magnitude of 
the tidal force (Equation |5]l. For the studies above, we chose 
?7f = 3, cf = 1, f = 9/ 2 and 7f = 2 for consistency 
with Villaver & Livio ' (2009). Now we explore the impact of 
changing these parameters to other values given in the litera- 
ture. 

The parameters governing the frequency dependence, cp 
and 7f, turn out to be relatively unimportant, because the con- 
vective time-scale is usually shorter than the tidal periods (at 
most tconv is 1.15 years, for the 5Mq star: see Figure [T2]|. 
Varying 7f while keeping the other parameters unchanged 
thus has no effect on our results. There is a small effect 
in changing cp: where we have picked cf = 1, comparing 
a complete tidal forcing period to the convective time-scale, 
other suggestions are for c f = 2 (Zahn 1966) or cp = 2n 
jGoldreich & Keeievlll977 ). Picking the larger as having the 
strongest effect, we find a decrease in the radius of the moritu- 
rus ultimus of around 0. 1 AU for a circular Jupiter and 0.2 AU 
for an eccentric Jupiter orbiting a 5 Mq star. The eccen- 
tric planet is more strongly affected since the eccentricity 
brings higher frequency Fourier components into pla y. This 
will make an application of the full theory of Eggleto n et alJ 
(tl998) for higher eccentricities impossible, as this theory is 
inconsistent with the responses to the Fourier components be- 
ing independent. Dealing with higher eccentricity planets will 
require explicit high-order eccentricity expansions of the tidal 
potential. 

There is a similar effect in changing the estimated convec- 
tive time-scale tconv with the parameter r?F. Changing r]p 
from our 3 to an alternative 1 (e.g., Zahn 1977) increases the 
convective time-scale by almost 50%, weakening tidal forces 
by increasing the denominators of equations ^ and (|2]l. In 
this case, the radius of the Jovian moriturus ultimus decreases 
by 0.1 — 0.2 AU, depending on stellar mass. We also reduced 
the multiplicative factor /' b y a factor 10, as suggested by 
the numerical simulations of Penev et al.l (l2009h . although we 
caution that these were for Main Sequence Solar-type stars. 
This resulted in a radical decrease in the radius of the Jovian 
moriturus ultimus, by 0.67 AU for the 1 Mq star and 1.25 AU 
for the 5 Mq star 

Finally, we note that the planetary tide had no noticeable ef- 
fect on the orbital evolution, even when the quality factor was 
very low: varying Qpj from 10^ to lO^'' for giant planets and 

from 10" to 10"^ for terrestrials had no effect on the outcome 
of the integrations. Planetary tides are unimportant due to the 
extremely small value of Rpi/a. 

5. DISCUSSION 

5.1. Relation to previous work 

Our study natura lly follows on from that of 
IViUaver & L ivio (2009), studying the fate of planets around 
RGB stars, i Villaver & Livio (2009 ) found that Jovian planets 
initially within 3 AU of a 1 Mq star would be drawn into the 
stellar envelope and engulfed. While this is immaterial for 
determining the maximum initial radius of engulfed planets 
on the AGB, it may affect the minimum orbital radius after 
the AGB, since the minimum final radius is attained for some 
minimum initial radius. If planets at this radius are engulfed 
on the RGB, the minimum post-AGB radi us may be larger. 
We no te, however, that in the integrations of IVillaver & Livid 
(|2009), planets were sometimes "stranded" by tidal forces, 
ending up inside the maximum RGB stellar radius due 
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Figure 12. Convective time-scale tconv over the thermally pulsing AGB life- 
time of a 5 M0 star. The parameter rjp = 3. 

to the small amounts of mass loss and the smooth stellar 
radius evolution. Furthermore, observations of Horizontal 
Branch stars have revealed some close companio ns that have 
appar ently survived t he RGB phase (ISilvotti et al.l 120071 : 
,Charpinetet all 120111: ISetiawan etaP 120101) . showing that 
some planets m ay be able to s urvive engulfment in the stellar 
envelope (e.g., ISoke3 [19981 iBear & SokeJ I20T2I) . Hence, 
the region closer to the star can be repopulated at the end of 
the RGB, and this makes our determination of the minimum 
post-AGB radius sound. For higher-mass stars this is not an 
issue, since any planets beyond a fr action of an AU are safe 
on the RGB (IVillaver & Livioll2009h . 

We also build naturally on the work of IVillaver & Liviol 
( 120071) . who studied planetary survival around AGB stars but 
neglected tidal evolution, on the grounds that the maximum 
stellar radius, and hence the strongest tidal force, exists only 
briefly at the AGB tip. It is clear from the above results how- 
ever, that for Jovian planets the tidal forces are indeed signif- 
icant over much of the AGB phase. IVillaver & Liviol (l2007h 
estimated the radii of the innermost surviving planets by as- 
suming that they would come from planets whose initial ra- 
dius was at the maximum stellar radius. However, since sig- 
nificant stellar mass loss occurs before the maximum radius is 
reached, the innermost survivor can originate from inside the 
maximum radius. Hence, the values for the minimum radii at 
which planets can be found on leaving the AGB phase in that 
work, ranging from 3 AU for a 1 Mq star to 30 AU for a 5 Mq 
star, are somewhat overesti mated: We find, on the contrary, 
2 and 10 AU respectively. IVillaver & Liviol (i2007) however 
treated another threat to planetary survival: atmospheric evap- 
oration due to the intense XUV radiation during the PN phase. 
They found that the gaseous envelopes of Jupiter-mass planets 
ending up within 3 — 5 AU will be totally evaporated during 
the PN phase, and hence that for lower-mass WDs that which 
determines the survival of Jovian planets is evaporation, not 
tidal evolution. Note however, that the cores of these planets 
will remain, and the lower-mas s planets we have considered 
will not be as strongly affected. IVillaver & Liviol (12007 ) also 
concluded that for higher-mass stars the survival of planets to 
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the WD phase is not affected by evaporation, due to their very 
large semi-major axes. Even with our significant reduction of 
the axis of the innermost survivor around these stars, we can 
agree that evaporation will not be significant for them. 

A study of the combined effects of tides and mass loss over 
the whole stellar life time has recently been carried out by 
iNordhaus et al.l(l2010l) . This study investigated the fates of Jo- 
vian planets, brown dwarfs and low-mass stellar companions 
to stars of mass 1 — 3 Mq. They compared the tidal formal- 
ism of Zahn and the Q formalism, finding that the latter yields 
considerably weaker tides. This may be the result of their us- 
ing a Q calibration for Main-Sequence stars, not giants, as we 
re marked in Secti o n[2l 

INordhaus et al.l (120101) also investigated a range of stellar 
mass-loss rates, varying the Reimers 77 parameter. However, 
the commonly-used Reimers (1975) rate does not apply to the 
AGB phase. As we have discussed in the previous sections, 
mass loss is a determining factor in the evolution of the orbit 
of the planet, and this is especially true on the AGB where 
large mass loss takes place. We therefore spend some time 
here discussing mass loss on the AGB and its effects on plan- 
etary orbital evolution. 

Renzini (1981) showed that a much stronger wind 
than the one allowed by the Reimers parameterization — a 
"superwind" — develops at a crucial point in the star's evolu- 
tion, reaching mass loss rates of 10^^ to 10^* M0 yr~^. The 
Reimers relation for the AGB phase does not agree with ob- 
servations through the entire AGB phase: it gives mass-loss 
rates that increase too slowly, take too l ong, and r each values 
too small (see e.g. Willson 2000; Bowen & Willson 199 iT) . 
Furthermore, Reimers mass-lo ss rates do not allow the forma- 
tion of planetary nebulae (e .g. lRenzini|[l981l : llben & Renzinil 
119831; ISchoenbernej|1983h . and fail to repro duce the high 
mass-loss rates observed in AGB stars (e.g. [Winters et al.l 
UOOO; Wood 1979), the increase in mass loss during the AGB 
ascent (Wood 1979; Willson 2000), and the observed C-star 
luminosity functions in the Magellanic Clouds and the ther- 
mally pulsing AGB l ifetimes of M and C stars in Magellanic 
Cloud star clusters (Mari go et ani2008l) . 

It is important to keep in mind that although mass loss is 
a crucial process in the evolution of stars along the AGB, 
it cannot be calculated from first principles since it requires 
the use of dynamical model atmospheres that consider time- 
dependent dynamics (shock waves and winds), radiation 
transfer (strong variable stellar radiation field), and dust and 
molecular formation processes. Several stellar evolutionary 
models that follow the temporal behaviour of the mass loss 
during the AGB together with th e thermal pulsation of the 
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Wachter et al.ii2002) . Although the mass-loss rates are not de- 



rived from first principles in these models, and most of them 
rel y either on the dyna mical mode l atmosphere calculations 
of iBowenI (il988 ) and lArndt et all ([T997i) . or onjhe semi- 
empirical mass-loss rate formula derivations of lWoodI (119901) . 
they do provide an opportunity to study the extensive history 
of mass loss on the AGB and beyond. The comparison of the 
different mass-loss prescriptions with observations of individ- 
ual stars in the AGB evolutionary phase is complicated given 
the variable nature of the star (and the mass loss), and a full 
discussion on the subject would be lengthy and beyond the 
scope of this work. However, the I Vassiliadis & WoodI (Il993h 
prescription is the most widely used, and not only has it not 
been ruled out by observations but it has even been favoured 
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Figure 13. Mass-loss rates for 5 Mq stars under our model (solid) and SSE 
models with r] = 0.5 (dotted), 1.0 (dashed) and 5.0 (dot-dashed). 



over other parame terizations (see, e.g., iMarigo et al.l 119981: 
'Zilistr a et all2002l) . 

No te also that the maxim um mass-loss rate adopted in the 
IVassiliadis & WoodI ([1993) models is limited to the radiation 
pressure limit. Dust-driven wind models allow values to be 
up to two times this limit. If the maximum mass-loss rates 
increase, the timescale of the evolution of the star is modified 
accordingly. By considering the evolution of the mass-loss 
evolution provided by a set of stellar evolution calculations, 
we are simulating a realistic scenario with no free parameters 
that allows a more self-consistent treatment of the problem. 

To compare our mass-loss prescription with the Reimers 
prescription, we show the mass-loss rate from our 5 Mq 
model in Figure [T3] We also show on this plot the mass-loss 
rates from models using the SSE code (Hurley. Pols. & Toi3 
2000) with Reimers mass-loss rates with parameter 77 — 0.5, 
1.0 and 5.0. Under the Reimers prescription, more mass loss 
occurs earlier on the AGB than in our chosen models regard- 
less of the chosen value of 77 in the SSE models. Since mass 
loss occurs earlier using the Reimers prescription, planets can 
move out of reach of tidal forces before these have a chance 
to act. 

In the context of comparing our results to those of 
INordhaus et al.l (!201 0|), we now note the following points. 
First, the INordhaus e t alj (120101) stellar models fail to repro- 
duce the initial-to-final mass relation using the Reimers pre- 
scription: note that their 3 M0 star ends the AGB with a 
mass of 1.1 or 1.0 Mq with 77 = 0.7 and 1 (J. Nordhaus, 
2012, private communication), with only 7? = 5 approach- 
ing the canonical value of O.TMq (e.g.. IWeidemannI [19871 ; 
iKalirai et al.l2008 V One should then expect that, using a com- 
parable formalism for the tides, ou r morituri ultimi should 
agree better with the Nordhaus et al.l (1201 0) determination of 
maximum semi-major axis that is tidally engulfed which uses 
the largest (and unrealistic) value of the mass-loss parame- 
ter 7y = 5. We find however the opposite; that is, our re- 
sults using the Zahn formalism give compa rable values for 
the radii of the morituri ultimi only for the INordhaus et alj 
(.2010) calculations with Reimers i] — 0.7 and 1, i.e., those 
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that give an unrealistic final mass. Second, our results using 
the Zahn formalism give somewhat larger values for the radii 
of the morituri ultimi for lower-mass stars (for Jovian planets 
around 1 Mq stars, 2.7 AU versus 1.6 AU) but similar values 
for higher-mass stars (for Jovian planets around 2.5 Mq stars, 
3.8 AU versus ^ 4AU). We believe that the reasons behind 
this large discrepancy are as follows. First, under the Reimers 
prescription, mass loss is continuous from early times, and 
hence planets' orbits can expand befor e tides have a chance 
to act, as described above. Second, the Nordhaus et alJ (120101) 
models have a higher final stellar mass and a smaller maxi- 
mum radius than ours. Both of these rn ean that tidal forces 
are weaker in the iNordhaus et alj ( 120101) study than in ours, 
and hence that the radius of the moriturus ultimus is lower. 
For higher-mass stars, however, the final stellar mass in the 
[Nordhaus et al. ( 2010) models is much larger than it ought 
to be, and therefore the orbital expansion due to mass loss is 
less. Hence both the tidal decay and the adiabatic orbital ex- 
pansion are weakened, which conspire to give a radius of the 
moriturus ultimus similar to ours. 

The stellar evolutionary model is one of two major sources 
of uncertainty in our calculations, the other being the tidal 
model. It may be asked which is the more important. For 
1 Mq stars, our moriturus ulti mus radius of 2.7 AU is consid- 
erably larger than those of N ordhaus et al.l (120 IQ) using the 
same tidal prescription (0.6 to 1.9 AU, 1.6 AU for the same fi- 
nal White Dwarf mass). The biggest single-parameter change 
in the tidal model — changing /' to align with the hydrody- 
namical simulations of Penev et al.l (12009) — gives a change 
of 0.7 AU for this stellar mass. Hence, it appears that, de- 
spite uncertainties in tidal theory, a significant improvement 
can indeed be obtained by improving the stellar models. 

5.2. Implications 

We have shown that, around WDs whose progenitor masses 
range from 1 to 5Mq, planetary companions may be ex- 
pected outwards of 1 .4 AU to 10 AU for Jovian planets around 
1 Mq Stars to 5Mq stars, and outwards of 2.5 AU to 15 
AU for Terrestrial planets around 1 Mq stars to 5 Mg stars. 
The closer planets may experience significant evaporation, 
yet their cores may remain. If engulfed in the stellar enve- 
lope, low-mass planets should not have much chance of sur- 
vival. However, planets of higher mass than those considered 
here, or their remains, might survive the Common Envelope 
phase and be found interio r to the limits provided in this work 
(e.g.. Bear & Sokeiil2012l) . The distant planets that we expect 
to survive may be de tec table by dire ct imaging, but only if 
massive enough (Burleig h et al.ll2002l) Tirning rnethods (e.g., 
Silvottietal .112011) and'transits (e.g.. lAgolll20lTt iFaedi etal] 
201 Ih are more sensitive to close-in planets, for which sur- 
vival of the AGB phase may be a challenge. However, we 
do not wish to discourage observers from trying to find plan- 
ets within these limits given that other mechanisms than those 
considered here might be at work and that successful discov- 
eries will only stimulate further theoretical work. 

Any planets that do survive the AGB phase can still 
have a significant dynamical influence on any other bod- 
ies such as remnant comets or asteroids in the system, and 
in this way can bring about the delivery of pollutant ele- 
ments to the neighbourhood of the WD. Such delivery mech- 
anisms require dynamical instability, w hether local instabil- 
ity of the metalliferous planetesimals dBonsor et al.l 1201 It 
iBonsor & W vatt"2012- Debes et al."2012') or wholesale global 
instability tDebes & Sigurdsson.200Z D. Veras et al., MN- 
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Figure 14. Change in location of Jovian planets (dashed lines) and secular 
resonances (sohd lines) around a 1 Mq star. 



RAS, submitted). 

It may be thought that the stability of the system that 
survives the AGB is enhanced by the effects of tidal 
forces destroying the inner planet, or separating it from 
those further out. Multipl e-planet scattering, as studied by 
iDebes & SigurdssonI (l2002h . acts on a time-s cale strongly de- 
termined by the separation of the planets CChambers et alj 
ll996tlFaber & Ouillenll2007l) . and the detachment or destruc- 
tion of an inner planet would stabilize a system liable to un- 
dergo scattering. And a planet scattering particles inwards 
from an exterior belt, as envisaged by Bons or et al.l (1201 ll) . 
can be drawn away from the planetesimals by the tides. 

However, in other cases the detachment of an inner planet 
may make a system more prone to instability. Multiple- 
planet systems with inclined binary companions can be stabi- 
lized by their mutual Laplace-Lagrange perturbations against 
the otherwise devast ating Kozai effect ( Innanen et al. ,19971 : 
ISaleh & Rasid 120091). a protectiv e mechanism possibly at 
work in 55 Cnc ('Ka ibet al.ll20Tlb . If an inner planet is de- 
tached, the coupling between the planets may weaken enough 
to trigger instability after the AGB phase. In any system, the 
detachment of an inner planet may cause divergent crossing of 
Mean Motion Resonances and the sweeping of secular reso- 
nances, both of which can pump up the eccentricities of plan- 
etesimals, and which are thought to have ha d an important 
effect on bodies in the early Solar System (e.g.. lTsiganis et al] 
2005; Brasser et al. 2009). Finally, the closer the inner planet 
to the WD, the easier can be the delivery of small bodies being 
scattered by chains of planets ( Bonsor & Wyatt 20T2h . 

As an example of the movement of secular resonances, we 
can consider a system with two planets of Jovian mass, one 
initially at 2.69 AU and one at 4AU, around a 1 Mq star. 
Their final semi-major axes are 1 .78 AU and 6.96 AU. As the 
system evolves, the secular resonances move from 5. 1 AU and 
7.2 AU to 13.9 AU and 14.0 AU (see FigurelHl. A final semi- 
major axis of 14 AU corresponds to an initial semi-major axis 
of 8 AU, and so the secular resonances move out faster than 
orbits expand due to adiabatic mass loss. Hence, any plan- 
etesimals initially between 5 and 8 AU will be swept by one 
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or both resonances, pumping up their eccentricity and poten- 
tially rendering them vulnerable to scattering at the end of the 
AGB itself or beyond. The full working out of the system- 
wide effects of the tidal evolution of inner planets must, how- 
ever, be deferred for future work. 

6. CONCLUSIONS 

We have studied the orbital evolution of planets around 
Asymptotic Giant Branch stars ranging from 1 M0 to 5 Mq as 
the star expands and loses mass, with detailed stellar models 
giving the stellar mass loss and radius evolution. We consid- 
ered planets of Jovian mass and lower, which have insufficient 
energy to unbind the stellar envelope and will likely be de- 
stroyed upon ent ering it. To compute the tidal forces we used 
the formalism of lZahnl([T977i) . appropriate for AGB stars with 
massive convective envelopes. Planetary orbital evolution is 
then a contest between tidal forces pulling the planet in and 
mass loss causing orbital expansion. 

Jovian planets experience strong tidal forces, and Jovian 
planets that begin on orbits outside the maximum stellar ra- 
dius can be drawn into the envelope. The maximum initial 
radius at which a Jovian planet on a circular orbit can be 
drawn into the envelope ranges from 2.6 AU for a 1 M0 star 
to 5 AU for a 5 M0 star Lower mass planets feel tidal forces 
less strongly, and planets that begin on orbits interior to the 
maximum stellar radius can escape engulfment due to their 
orbits expanding under stellar mass loss: Terrestrial planets 
on initially circular orbits can be drawn into the envelope at 
initial radii of up to 1.5 AU to 2.8 AU over our mass range. 

The inclusion of a planetary eccentricity increases the max- 
imum initial radius for which planets can enter the stellar en- 
velope. In the case of Jovian planets this is due to the strength- 
ened tidal forces: the lifetime of an eccentric planet is reduced 
by about 10-20% relative to a circular one, and eccentricity 
decays almost to zero as the planet plunges into the envelope. 
Eccentricity is damped entirely by the stellar tide and plan- 
etary tides are ineffective. The eccentricities of Terrestrial 
planets barely decay, and their having smaller pericentres rel- 
ative to circular planets means that the maximum initial radius 
at which planets will be engulfed is much larger, and their life- 
time much shorter, by as much as one half. 

We then found the minimum semi-major axis at which plan- 
ets may be found around White Dwarfs after the AGB stage 
has ended. This is larger for more massive stars, both be- 
cause larger stars lose a greater fraction of their mass on the 
AGB, causing more adiabatic orbital expansion, and because 
for the larger stars the initial radii of the survivors are larger 
Assuming that any planet of Jovian mass or less entering the 
stellar envelope is destroyed, we do not expect to find Terres- 
trial planets within 2 AU of WDs of progenitor mass 1 Mq or 
within 15 AU of WDs of progenitor mass 5 Mq. Jovian plan- 
ets, which can have their adiabatic orbital expansion retarded 
by tidal forces even as they avoid the stellar envelope, may 
be found somewhat closer, at 1.5 AU or lOAU respectively. 
However, a certain fine-tuning is required in order that such 
planets avoid plunging into the envelope while still experi- 
encing a non-negligible tidal force: there is a region around 
0.5 AU in width outside the orbit of the maximum orbit that 
gets engulfed in which Jovian planets still experience strong 
tidal orbital decay. 

We studied the effects of changing the several parameters in 
the Zahn tidal theory. The systems we studied have convective 
time-scales a few times shorter than orbital time-scales, and so 
the unknown nature of the frequency dependence of the tidal 



forces are relatively unimportant. However, uncertainties in 
the overall magnitude of the tidal response can have a signif- 
icant effect, with a reduction of a factor 1 to align th e Zahn 
theory with the numerical simulations of iPenev et alJ (120091) 
reducing the radius at which Jovian planets can be engulfed 
by around 1 AU. 

Finally, we noted that in multi-planet systems the differen- 
tial orbital evolution of more and less distant planets could 
have significant implications for the stability of the system on 
the White Dwarf phase, and therefore for the pollution of the 
White Dwarf's atmosphere and its environment by destroyed 
planets or planetesimals. 
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